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I.  INTRODUCTION 


The  flow  of  liquid-vapor  mixtures  occurs  in  many  engineering 
applications,  such  as  drains,  the  handling  of  refrigerants  and  condensed 
gases,  blow-offs  from  turbines  or  boilers,  spacecraft  propulsion  and 
attituae  control  systems,  and  in  the  coolant  system  design  for  nuclear 
power  plants. 

I'o  accurately  analyze  a two-phase  flow  it  would  be  advantageous  to 
know,  a priori,  whether  the  flow  is  critical  or  subcritical.  The 
ability  to  predict  both  critical  and  subcritical  flow  rates,  as  well  as 
other  dynamic  and  thermodynamic  properties  of  the  flow  field,  is  impor- 
tant to  the  understanding  of  these  problems. 

Two-phase  flow  has  been  treated  extensively  in  the  literature,  due 
primarily  to  the  advances  in  nuclear  reactor  technology.  The  emphasis 
of  the  two-phase  flow  research,  however,  has  been  to  predict  the 
critical  flow  rate  and  the  pressure  wave  projKigat  i on . (!omprcssible, 
subcritical  two-phase  flow  has  received  relatively  little  attention, 
probably  because  maximum  or  critical  discharge  rates  resulting  from 
ruptured  steam  lines  are  of  primary  importance  in  the  design  of  nuclear 
power  plant  coolant  systems.  Subsonic  flow  is  usually  treated  as 
incompressible.  The  Bibliography  bears  witness  to  these  facts. 

Two-phase  flow  analyses  arc  complicated  by  dcjiartures  from  equili- 
brium (metastability),  slip  between  the  two  jihases  (non-liomogeneity)  , 
and  an  undefined  sonic  or  critical  flow  velocity  concept.  Also,  two- 
phase  flow  analyses  usually  depend  on  the  quality  range,  void  fraction 
range,  and  the  flow  regimes  under  investigation.  Mul t i -component , 
two-phase  flow  analyses  introduce  additional  complications.  A complete 
description  of  the  two-phase  flow  phenomenon  should  include  an  accurate 
subsonic,  compressible  flow  model;  a critical  jiressure  ratio  or  other 
criterion  to  describe  the  onset  of  critical  flow;  a critical,  compres- 
sible flow  model;  and  a jiressure  wave  jirojiagat  ion  jirediction  comjiatible 
with  the  subsonic  and  critical  flow  models.  These  models  should  be 
dejiendent  on  the  flow  regimes  defined  by  the  void  fractions,  sliji 
ratios,  mass  flow  rates,  and  (jualities  associated  with  these  flows. 

The  existing  two-phase  flow  models  are  primarily  designed  to  jiredict 
the  critical  flow  rate  of  a steam-water  mixture  discharging  from  a 
nozzle,  short  pijie,  or  orifice  [1-27].  Tlie  sonic  velocity  or  jiressure 
wave  propagation  in  a two-phase  flow  has  been  studied  in  |,i,  27-48,  14- 
17,  21,  2.3,  24).  Items  (7,  4(i,  47)  are  recent  survey  articles  sujijilc- 
mented  by  the  Bibl  iograjihy  given  here. 
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To  date,  the  models  proposed  by  Fauske  or  Henry  [5,  7,  8,  9,  11, 

12,  18,  25,  26,  31,  49]  or  those  proposed  by  Moody  [1-3]  are  the  most 
popular  methods  for  predicting  two-phase,  one-component  (steam-water) 
critical  flow  rates.  The  Aisch  model  [22]  is  capable  of  predicting 
either  critical  or  subcritical  flow,  and  liquid  effects  are  included 
through  an  isentropic  specific  heat  ratio  and  in  the  density.  Large 
scale  computer  codes,  such  as  RELAP4  [50]  calculate  either  subcritical 
or  critical  flow,  but  this  is  accomplished  by  taking  the  lesser  of  several 
calculated  flow  rates,  for  a conservative  (low)  estimate  of  the  flow  rate. 
Thus,  the  onset  of  critical  flow  is  not  predicted. 


II.  THE  VENT  FLOW  MODEL 

A two-phase,  single-component,  homogeneous,  frozen,  adiabatic  vent 
flow  model  is  presented.  This  model  is  applicable  to  problems  associated 
with  the  postulated  rupture  of  high  energy  pipe  lines  in  nuclear  power 
plant  coolant  systems  and  the  resulting  subcompartment  pressure  buildups. 
Items  [51  and  52]  of  the  Bibliography  provide  a detailed  description 
of  the  overall  problem.  This  analysis  considers,  in  effect,  only  one 
component.  Multi -component  (several  vapors  or  gases)  can  be  incorporated 
into  the  analysis  by  appropriately  weighting  the  gas  components  by 
standard  methods  such  as  discussed  by  [22]. 

At  any  given  time  during  the  transient,  the  total  conditions  in 
the  nodes  upstream  and  downstream  of  the  vent  are  assumed  given,  see 
[52],  for  example.  Then,  the  subcritical  or  critical  flow  that  will 
occur  across  the  vent,  with  a specified  area  and  resistance  coefficient, 
as  a result  of  the  imposed  pressure  differential  can  be  determined.  This 
mass  flow  rate  is  then  used  to  update  the  node  inventories,  and  the 
procedure  is  repeated. 

The  vent  flow  model  considers  two  cases.  First,  isentropic  inlet 
effects  are  included  in  the  flow  model  when  the  vent  is  best  represented 
by  an  area  reduction  (contraction).  Second,  when  a contraction  of  the 
flow  does  not  occur,  such  as  in  subcompartment  nodal ization  studies, 
the  isentropic  inlet  effects  are  not  included. 
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At  the  onset,  the  following  assumptions  are  made; 

1.  The  flow  is  quasi-steady,  i.e.,  the  flow  at  any  point  in  time  is 
assumed  steady. 

2.  The  flow  is  one-dimensional. 

3.  The  flow  is  homogeneous. 

4.  The  flow  is  adiabatic. 

5.  No  mass  transfer  occurs  between  phases, 

6.  The  vapor  phase  of  the  mixture  is  a thermally  and  calorically 
perfect  gas. 

Pressure  changes  within  the  vent  due  to  gravity  are  negligible. 

8.  The  vent  loss  coefficient  is  constant.  However,  a two-phase  flow 
multiplier  could  be  employed,  or  any  two-phase  friction  coefficient 
model  which  can  be  represented  in  terms  of  the  unknown  flow  field 
parameters  could  be  used.  Thus,  this  assumption  may  be  removed  if  a two- 
phase  flow  loss  coefficient  expression  is  known.  Also,  vena  contracta 
losses  are  included  in  the  vent  loss  coefficient. 

9.  The  vent  area  is  constant  during  a time  step. 

10.  The  liquid  phase  is  incompressible  over  the  specified  range  of 
pressures  and  qualities. 

A geometric  representation  of  the  flow  field  under  discussion  is 
depicted  in  Figure  1. 

The  continuity,  momentum,  energy  and  state  equations  are  solved  for 
a two-phase,  single-component  mixture  using  the  assumptions  given  earlier. 
Algebraic  equations  are  generated  which  allow  the  determination  of  the 
vent  inlet  and  outlet  Mach  numbers  based  on  knowledge  of  the  total 
pressure  ratio  across  the  duct  and  the  vent  friction  factor.  Knowledge 
of  the  inlet  and  exit  Mach  numbers  in  turn  allows  the  calculation  of 
other  flow  field  parameters.  This  model  also  enables  critical  flows  to 
be  calculated  based  on  the  criterion  that  critical  flow  coincides  with 
sonic  flow  at  the  duct  exit. 


The  governing  equations  follow. 
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Figure  1.  Two-Phase  Flow  Through  A Constant  Area  Duct 
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The  stagnation  enthalpy  is  conserved,  i.e., 


•Ih 


du^ 
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= dh  = 0, 
o 


H) 


where 

dh  = Xdh  + (1-X)  dh,-. 

8 i 


The  vapor  phase  is  an  ideal  gas. 


so  dh  = C d'P  . 
g P g 


The  total  derivative  of  the  liquid  phase  enthal])y  is 


because  the  liquid  is  incompressible  and  adiabatic  flow  is  assumed. 

The  flow  is  assumed  homogeneous,  so  the  energy  equation  (IJ 
becomes 


XC  dT  + = 0, 

P,  g 2 


(2) 


where  the  quality,  X,  is  defined  by 


m 

X = 


m +mj. 
g t 


The  specific  volume  is  given  by 


or 


V = Xv  + (l-X)v,- 
g ^ 


for  Vg  >>  v^.  In  fact,  this  approximation  is  acceptable  (to 
within  5%)  for  qualities  above  0.2  and  pressures  below  2070  k Pa. 
The  equation  of  state  for  an  ideal  gas  is 


Pv  = RT 
g 


or  Pv  = XRT.  (4J 

The  sonic  velocity  of  a homogeneous,  frozen  two-phase  mixture, 
neglecting  the  compressibility  of  the  liquid  is  [..ij 


a (Pf/Pg) 


Introducing  the  void  fraction,  a,  of  a homogeneous  two-phase  flow 
mixture , 


(l-X)p  +Xpj, 


the  sonic  velocity  becomes. 


a^p  = X a^2  for  p^  « p^. 


The  two-phase  sonic  velocity  for  an  ideal  gas  becomes 


a2p  = XyRTg  = yPv. 


/ 


'4 


Thus,  the  homogeneous  two-phase  Mach  number  may  be  expressed  as 


m2 


XyRT 


f9) 


Thus,  equation  (2)  becomes,  for 


C 


P 


YR  > 
Y-1 


dT 

T 


+ 


0. 


flO) 


For  the  flow  through  a duct,  the  momentum  equation  may  be  written 


oudu  . dp  • = 0,  no 

where  f is  the  resistance  coefficient  (or  four  times  the  friction 
coefficient  defined  by  Shapiro  [53]). 

The  continuity  equation  is 


C.  = pu, 


or 


1 du2  dp 


0. 


i:i 


Using  the  definition  of  Mach  number,  (9),  the  momentum  qu.it  ion 
(11)  may  be  rewritten  as 


yM2  du2  ^ d]l  + 1 V,?  liL  = (( . ( I 'M 

2 y r 2 1' 


The  equation  of  state,  (4)  may  be  expressed  as 


+ ill  n ») 

P ' p T 
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and  (10)  becomes 


^ ^ ^ Yj_  m2 


du^ 


= 0. 


Equations  (12)  and  (15)  may  be  combined  to  yield 

- dlVl’ 


du^ 

..2 


1/2[1+(y-1)M2) 


Thus,  (13)  may  be  written  as 

dP  _ -yM^[U(y-1)M^J  fdL 


2[1-M2] 


[) 


(15) 


(16) 


(17) 


Now  returning  to  (9)  and  using  logarithmic  differentiation, 

(18) 


dM^ 

m2 


du2 

..2 


£T 

T 


or 


dM^ 


= ^ . I-  "1  • 


using  (10) 


du2 

m2  u‘ 

Equations  (18)  and  (16)  may  be  combined  to  yield 


P 


1+(Y-1 JM2 
2+(y-1)M2 


dM2 

m2 


and  equations  (19)  and  (17)  yield 

,2  [ 1 * ^ m2 


dM2 

m2 


yM^ 


1 - m2 


fdL 

[) 


(19) 


(20) 
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The  equations  derived  above  are  identical  in  form  to  those  obtained 
for  single  phase  flow  through  a constant  area  duct  f53j.  However,  the 
form  is  only  cosmetic,  since  the  Mach  number  has  been  redefined  according 
to  (9)  to  account  for  liquid  effects. 

Equation  (20)  may  be  integrated  between  the  duct  inlet  (station  1) 
and  exit  (station  2),  to  yield 


fL 

D 


Also,  (19) 


m2  - m2 
2 1 

+ m 1 

(1* 

^ M^) 
2 2 

m2 

1 

y 

m2  m2 
. 1 2 ^ 

' -y  X 11 

2y 

(1  + 

2 1 

m2 

2 

integrated  between  the  same  two  stations  yields 


(21) 


(2+(y-1)M2)m2 
1 1 


(2+(y-1)m2)M? 

9 ? 


1/2 


(22) 


Critical  flow  occurs  when  the  vent  exit  Mach  number,  M^,  equals 


one . 


The  mass  flux  per  unit  area,  G,  may  be  written  as 


G = P u = M /yf’  P . (23) 

11  111 


using  (9)  and  the  continuity  equation. 

The  above  expression  for  G is  applicable  for  the  flow  between  two 

reservoirs  where  an  area  reduction  (contraction)  is  not  present,  i.e., 

wall  shear  represents  the  only  frictional  loss.  In  this  case,  the  flow 

from  station  1 to  station  2 is  accelerated  by  the  pressure  difference 

P , > P ... 
ol  exit 
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When  an  area  reduction  occurs  between  the  inlet  reservoir  and  the 
vent,  as  shown  in  Figure  1,  an  isentropic  inlet  effect  is  included. 

In  this  case,  (2)  integrated  between  the  reservoir  and  vent  inlet  yields 


XC  (T  - T ) = u2/2  (24) 

P 01  1 1 


where  u^j  ^ 0,  for  the  static  reservoir. 


Since  the  flow  is  adiabatic  and  homogeneous,  the  gaseous  phase 
will  expand  isentropical ly  where 


or 


Pv  ^ = constant, 

g 

Y 

Pv  = constant, 


(25) 


from  (3)  for  frozen  flow.  Of  course,  these  constants  are  not  the  same. 
Equations  (4)  and  (25)  yield. 


and  thus,  (24)  becomes 


(26) 


T 

1 


- 1 


Y-1 

2 


yxrt 

1 


or 


P 


01 


P 


1 


1 


+ 


7-1 

2 


m2 


(27) 
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The  back  pressure  in  the  sink  node,  P , , is  assumed  to  be  the  exit 
' ’ exit 

static  pressure,  P , where  the  sink  node  represents  a plenum. 

? 

The  pressure  ratio  (22)  then  becomes. 


ex  1 1 


ol 


(2+(y-1)M2)m2 
1 1 

(2+(y-1)M2)m2 
2 2 -■ 


H 1/2 


(1 


2 1 


-Y 

Y-1 


(28) 


For  isentropic  flow  between  the  reservoir  and  vent  inlet,  using 
the  isentropic  flow  relationships. 


G = 


(l*  m2  ) Y-1 


^ Pois’d! 


1 + 


Y-1 


m2 


1/2 


(29) 


Also,  several  auxiliary  equations  may  be  derived  to  provide 
additional  vent  flow  data. 


Using  the  expression  for  dP/P  from  equation  (17),  equation  (!(>] 
becomes 


du  _ yM^  fdl. 

2(1 -M^) 


( .^0) 


1 V 
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-i =- 


■A 


•d 


and  (17)  and  (19)  yield 


dM^  _ [2-i-(y-1)M^]yM^  fdL 

m2  2(1-M2)  D 


(31) 


Also,  equations  (30),  (31)  and  (18)  imply 

fdL  . 
> 

n 


^ 2(1 -m2) 


(32) 


and  from  <50)  and  (12), 


dp  _ yM^  fdL 
^ 2(1-M2)  1) 


Next,  equation  (27)  implies 


(33) 


dP 

o 

P 

o 


dP  yM^  dM2 

’’  2fl+  ^ m2]  m2 


or  from  (17)  and  (31), 


o yM^'  fd  L 
P ^ ■ 2 1)  ■ 


(34) 


(35) 


All  of  these  equations  reduce  to  the  conventional  gas  dynamics 
equations  [53]  for  single-phase  compressible  flow,  when  the  quality  is 
equal  to  one. 

Using  (31)  to  relate  the  friction  factor  to  the  Mach  number,  (30), 
(32),  (33)  and  (34)  can  be  integrated  to  yield: 
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u 

2 

u 

1 


(2+Cy-1)M2 
2+(y-1)M2 


(36) 


T 

2 

T" 

] 


2 


2+(y-1)M 

1 


2+(y-1)M2 

2 


(37) 


and 


M 

1 

nT 

2 


-I 

2+(y-1)M2  2(y-1) 

2 

2+(y-1)M2 

1 


(58) 


Equation  (38)  is  included  only  for  the  sake  of  completeness;  (28)  is 
used  for  the  total  pressure  ratio  for  the  actual  vent  flow. 

The  sonic  conditions  are  defined  to  occur  when  the  exit  Mach 
number  reaches  one.  The  critical  pressure  ratio  becomes 


When  the  given  pressure  ratio,  less  than  the  critical  pressure 

ratio,  the  flow  is  taken  to  be  critical. 


III.  SOLUTION  TECHNIQUE  AND  RESULTS 


Results  are  obtaned  by  the  simultaneous  solution  of  equations  (21) 
and  (22)  or  (28)  for  and  M^.  Equation  (22)  is  used  when  inlet  effects 

are  ignored,  e.g.,  there  is  no  connecting  duct  between  stations  one  and 
two.  In  this  case,  the  stagnation  pressures  are  used  in  place  of  the 
static  pressures  in  equation  (22),  and  the  flow  is  artificially  accelerated 
between  stations  one  and  two.  Equation  (28)  is  used  when  inlet  effects 
are  included,  or  a flow  restriction  (or  duct)  exists  between  stations  one 
and  two.  The  total  pressure  in  the  nodes  upstream  and  downstream  of 

the  vent  is  assumed  known.  The  friction  factor  is  given  or  can  be 

expressed  as  a function  of  the  flow  field  parameters.  This  functional 
form  of  the  friction  factor  can  then  be  expressed  in  terms  of  the  inlet 
and  exit  Mach  numbers  using  the  relationships  given  previously.  The 
governing  equations  are  then  more  complicated,  but  can  still  be  solved 
simultaneously  for  and  M2.  The  two  governing  equations  are  reduced 

to  one  nonlinear  algebraic  equation.  The  critical  flow  equations  are 

solved  for  M2  equal  one  and  the  resulting  critical  pressure  ratio  is 

calculated.  This  critical  pressure  ratio  is  compared  to  the  actual 
pressure  ratio  to  determine  if  the  flow  is  choked.  If  the  flow  is 
choked,  the  mass  flux  is  calculated.  If  the  flow  is  subsonic,  the 
governing  equations  are  solved  simultaneously  for  M^  and  M2,  and  then 

the  mass  flux  is  calculated.  Once  the  inlet  and  exit  Mach  numbers  are 
known,  other  flow  field  parameters  can  be  readily  calculated. 

The  actual  solution  is  achieved  by  attempting  to  solve  the  final 
nonlinear  equations  using  the  Newton- Raphson  method.  If  this  method 
does  not  converge,  solve  the  equations  by  the  Bisection  method. 

Several  parametric  studies  utilizing  this  flow  model  have  been 
performed  and  inconsistent  behavior  was  not  observed.  Figures  2A,  2B, 

2C  and  2D  present  the  mass  flow  rate  versus  inlet  pressure  for  a 

specified  back  pressure.  The  quality  and  friction  factor  were  varied. 

A steam-water  mixture  was  considered  with  y=1.1,  although,  in  general, 

Y is  calculated  from  the  known  reservior  data.  Figures  2A  and  2B  include 
the  isentropic  inlet  effect;  Figures  2C  and  2D  do  not. 

Figures  3A,  3B  and  3C  compare  the  critical  flow  rates  versus 
critical  pressure  to  the  analytical  data  of  Fauske  [8]  for  steam-water 

mixtures.  Two  of  Fauske's  models  are  shown  [8J.  Comparison  is  also 

made  to  the  experimental  data  of  Faletti,  Moy,  Fauske  [8],  also  [4), 

[10];  and  Klingebiel  [13].  Figures  3 indicate  good  agreement  with  the 
experimental  data  and  the  theoretical  predictions  of  Fauske,  especially 
at  high  qualities  (X>.5).  The  results  of  this  model  also  agree  very 
well  with  Moody's  theoretical  results  as  presented  in  [47].  For  a 
quality  of  one,  the  equations  associated  with  the  vent  flow  model  reduce 
to  the  classical  gas  dynamics  expressions  [53],  as  expected.  Thus, 
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the  theory  is  exact  for  a quality  of  one;  excellent  for  qualities  above 
0.5;  acceptable  for  qualities  as  low  as  0.2;  and  poor  for  qualities  less 
than  0.1. 


Finally,  Figure  4 illustrates  the  variation  of  vent  inlet  pressure 
with  mass  flow  rate  over  a range  of  qualities. 

The  results  of  the  vent  flow  model,  including  isentropic  inlet 
effects  exhibit  excellent  agreement  with  Aisch  [22J  for  frictionless 
flow  using  a constant  y(=1.1).  This  model  differs  from  that  of  Aisch  in 
the  treatment  of  the  friction  coefficient  (it  is  an  inherent  in  the 
equations  of  the  vent  flow  model;  Aisch  modifies  the  calculated  flow 
rate  by  use  of  a discharge  coefficient]  and  in  the  inclusion  of  liquid 
effects  (Aisch  considers  an  isentropic  flow  and  modifies  the  specific  heat 
ratio  and  the  density  to  account  for  two-phase  effects;  in  the  vent  flow 
model,  frozen  flow  is  assumed  and  two-phase  flow  effects  are  included  in 
the  sonic  velocity,  density,  and  specific  heats). 

IV.  CONCLUSIONS 

A two-phase  flow  model  is  presented  which  predicts  the  critical 
pressure  ratio  for  choked  flow,  calculates  both  subsonic  and  critical 
flow  rates,  calculates  other  flow  field  parameters  of  interest,  and 
utilizes  a sonic  velocity  equation  consistent  with  the  flow  field 
model.  The  model  can  he  readily  extended  to  include  mult i -component 
two-phase  flow  and  a flow  loss  coefficient  dependent  on  the  flow  field 
parameters  (or  including  a two-phase  multiplier).  The  model  is  limited 
by  the  assumptions  of  no  momentum  exchange  between  phases  and  no  mass 
or  heat  transport,  i.e.,  a frozen,  adiabatic,  homogeneous  flow  is 
assumed.  However,  these  assumptions  tend  to  offset  one  another  and 
the  predicted  results  show  good  agreement  with  the  existing  critical 
flow  data.  Ihe  vent  flow  model  is  cajiable  of  providing  useful  design 
analyses  in  a high  quality  (X>..t),  high  void  fraction  (i>.9)  flow- 
regime  . 
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a Sonic  Velocity 

a Void  Fraction 

Cp  Coefficient  of  Specific  Heat  at  Constant  Pressure 

dL  Increment  of  Vent  Length 

D Vent  Hydraulic  Diameter 

f Resistance  Coefficient 

fL  = K Friction  Factor 

iT 

G Mass  Flow  Rate 

h Lnthalpy 

y Specific  Heat  Ratio 

L Vent  Length 

M Mach  Number 

m Mc'i  s s 

P I’ res  sure 

R Gas  Constant 

T Temperature 

p Density 

u Velocity 

V Specific  Volume 

X Quality 
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LIST  or  SYMBOLS  (Cont . ) 


Subscripts 
g Vapor 

f Liquid 

1 Vent  inlet 

2 Vent  exit 

0 Stagnation  or  Reservoir 

TP  Two-Phase 

Superscripts 

* Sonic  or  Critical  Condition 
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1 Commander 

US  Army  Missile  Research 
and  Development  Command 
A'lTN:  DRDMI-R 

Redstone  Arsenal,  AL  35809 

1 Commander 

US  Army  Tank  Automotive 
Development  Command 
ATTN ; DRDTA- RWL 
Warren,  MI  48090 


No.  of 

Copies  0 r ^ ao  i^aj:  ion 

2 (T)mmander 

US  Army  Mobility  licpiipment 
Research  t,  Development  (imd 
ATTN:  lech  Docu  Cen , Bld}>  31 

DRSMi;-RZT 

fort  Belvoir,  VA  22060 
1 Commander 

US  Army  Armament  Materiel 
Readiness  Command 
ATTN:  DRSAR-Llif 

Rock  Island,  IL  61202 

1 Commander 

US  Army  Armament  Research 
and  Development  Command 
ATTN;  DRDAR-LC,  Dr . .1 . fras  i or 
Dover,  N.J  07801 

1 Commander 

US  Army  Harry  Diamond  Labs 
A'lTN:  DRXDO-Tl 

2800  I’owder  Mill  Road 
Adel  phi,  MD  20783 

1 Director 

US  Army  TRADOt,'  Systems 
Analysis  Activity 
ATTN:  ’ AI'AA-SA 
White  Sands  Missile  Rans;e 
NM  88002 

1 University  of  Illinois 
Department  of  Aeronautical 
and  Astronaut!  cal  lin^  . 

ATTN;  Professor  S.  M.  Yen 
Ur ban a,  IL  61801 


Aberdeen  Crmui^ 

Marine  Corps  l.n  Ofc 
Dir,  USAMSAA 
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